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STATISTICAL MECHANICS OF SUPERNOVAE 


Frank C. Jones 
ABSTRACT 

In most treatments of the supernova origin theory of cosmic rays, average 
values of such observables such as cosmic ray flux, bulk streaming velocity, 
total energy, etc., are calculated (or estimated) and compared with observation 
to determine the parameters of the theory. Since supernova explosions are es- 
sentially statistical events the question arises as to whether or not the observed 
values have any simple relationship to average values; in other words are fluc- 
tuations important? To investigate this question I have considered that the 
variables describing the supernova injection event (position, time, etc.) are random 
variables and that our galaxy is a sample from an ensemble of galaxies. In this 
manner one can, in principle, calculate all higher moments of the cosmic ray 
observables. It turns out that an important parameter may be derived for any 
particular model of cosmic ray transport. If a single supernova event fills an 
effective volume V eff with cosmic rays for a time r c and there are n super- 
nova events per unit time per unit volume on the average, the quantity (n) 
nT c V eff becomes essentially an inverse "discreteness parameter" in the sense 
that when it is large fluctuations are relatively unimportant, but when it is small 
fluctuations can dominate the situation. For currently considered models of our 
own galaxy the situation appears to be borderline. 
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I. Introduction 


For some time now supernova explosions have been considered a leading 
contender for the role of primary source of galactic cosmic rays. This is based 
primarily on the fact that their rate of energy release appears to be sufficient 
to meet the requirements of the cosmic ray source 1-6 10 41 ergs/sec) and 
more recently on the possibility that a detailed model of the supernova explosion 
might predict the energy spectrum of cosmic rays. 6-8 

A second major theme in recent cosmic-ray research has been the question 
of how cosmic rays propagate through the galaxy once they have been produced 
by their source. This question has been approached from two different aspects; 
first what combination of source characteristics and average path length in the 
galaxy will produce the abundances of the nuclear species observed in the 
cosmic- rays. There is, of course, a vast literature on this question and I cite 
here only a few of the most recent articles. The second approach to this ques- 
tion is concerned with the fact that the bulk streaming velocity of cosmic rays in 
the vicinity of the earth appears to be very small (-€500 km/sec). A consid- 
erable discussion has appeared in the literature concerning how the structure of 
the galactic magnetic field 14 * 15 » 16 and various possible types of plasma instabil- 
ities 17 » 18 » 19 could be responsible for this fact. A common feature of the 
theoretical work that has been done so far in these matters is that for any of the 
various models considered average values of observable quantities are calculated 
and then compared with observation, in so far as possible. It is the purpose of 
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this paper to show that if cosmic rays do originate in supernova explosions the 
system we are observing is an essentially noisy one and the values of observable 
quantities that we measure may not have a great deal to do with the average 
values. In other words the probability that we are observing a rather large 
fluctuation away from the average may not be negligible. 

II. Statistical Method 

We begin bv considering the general response function for a supernova 
explosion. The density of cosmic rays at a point X and at time t due to the i th 
supernova explosion shall be designated by P^X, t). Since we shall consider 
the propagation characteristics of the galaxy to be uniform in time and one 
s 'atial direction we have p { (X. t ) = p(X - X . , t - t ( ) where p is some uni- 
versal function and X ; and t . are the time and position of the i th supernova 
explosion. The reason for considering only one spatial coordinate will be made 
clear in Section in. when we discuss the actual propragation model that we will 
use. 

In a like manner we may discuss the directional streaming velocity ot the 
particles from the i th event v(x~X i .t-t i ) s v 4 (X, t) where the directional 
flux of particles is given by j . (X. t ) = v. (X. t) p i (X, t). 
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At any point in space and time the total density of cosmic ray particles is 
given by 


P 



Pi 


where the sum is over all supernova events. For reasons of later mathematical 
simplicity we shall consider the sum to be over all past and future events and 
build all causality requirements into the response function p(X, t). In a similar 
fashion we have 


j 



i 


The bulk streaming velocity that would be observed is given by 



and since the age of all particles from the 
of all particles is just 




i th event is just t - t . the a x "irage 



Until now there has been no mention of statistical notions; the eoordi’iates 


of the supernova events, X if t . have been considered to have determined values 



and the various observable quantities to be exactly calculable, at least in princi- 
ple. In reality, of course, we have no idea of the values of the event coordinates 
(with possibly a few’ exceptions) that would be required to calculate the values of 
the observable quantities. Even if they were known to a reasonable degree of 
accuracy the subsequent computations would be extremely difficult if not 
impossible. 

This situation is reminiscent of the one encountered in statistical mechanics 
and, in fact, we shall treat this problem from a statistical mechanical point of 
view. We shall consider our galaxy to be just one sample of an arbitrarily 
large ensemble of galaxies. In this ensemble the event coordinates X ; , t . shall 
be considered to be random variables, distributed with some probability distri- 
bution P(Xj , t j , X 2 , 1 2 , ••• X { , t . •••). Our observable quantities p, j , v , 
etc. , are now functions of many random variables and while we will not calculate 
these functions themselves w r e may in principle calculate average values for any 
combination of them. 

We shall use the method employed by Rice^O in calculating electrical shot 
noise. It makes use of the idea of an ensemble average of any function of the 
random variables; If F is a function of the X it t . then we have 


<F> s 



P(Xi. *,) r(x,. *,) 


(i) 
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Since we are interested in - 00 < t i < + ® there are an infinite number of 
events and it is not clear what we want as our probability function P(x. , t . ). 
Before approaching this question we shall make some simplifying assumptions. 
First we shall assume that separate events are statistically independent, second 
we shall assume that the probability of an event occurring in an increment of X. 
is L' 1 dX f if X a < X. < X b where L = X b - X. and is zero otherwise, and 
third we assume that the probability of an event occurring in an increment of 
time dt . is /idt . where p. is a constant. 

We shall now consider the calculation of averages in three stages. We shall 
first consider the contribution to an observable quantity from only those events 
occurring in the time interval - T/2 < t . < T/2 where T is a large but 
arbitrary finite length of time and further we shall consider only that subset of 
the ensemble for which exactly n events take place in this interval.- For this 
subset the averaging operation has the form 

n /» X b -t/2 

< >„ T = L-» T-" n dX, dt, . (2) 

T‘ 1 Jx > J - T/2 

This operator is well defined and is concerned only with finite quantities. It 
may therefore be applied directly to any observable quantity in a straightforward 
manner. 

In the next stage we simply note that the quantity n will have a Poit son 
distribution over the entire ensemble so our average may be extended ic ».ne 
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full ensemble by suitably averaging over n. Our results will now be functions of 
the quantity (o)/T the average number of supernova events in the time interval 
T. Passing to the limit is achieved by simply noting that we have assumed 

to be a constant therefore independent of T. We simply replace it with the 
average number of events per unit time n. By this method we readily obtain 

/L) f dX, [ dt,A>, (X,, tj . 

J ~* 

P) 2 + (n/L) f dX i [ dt . Pj 2 (Xj , tj (3) 

etc. 

We shall now anticipate a result of the next section. In Section IB we will 
see that our response function p(x - X i , t - t .) is actually a function of the 
dimensionless parameters (x-xJ/L,. and (t-t.)/r c where L c and r c area 
characteristic length and time respectively. The integrals may therefore be 
written as 




etc. If we write (n) e nL c t c /l and (pA. (pj 2 ) . etc., as the integrals over 
the dimensionless parameters of p x , p 2 , etc., w'e may write equations (3) as 


<P> = <N> <P,> . <P 2 > = (N) 1 ^,) 2 * (NXP, 2 ) 




8 


etc. The parameter (n) is the average number of supernova events per unit 
characteristic length, per unit characteristic time, and since it is essentially 
the average number of events that are sensibly contributing to the state of 
cosmic-ray-affairs at any one point in space and time it plays the role of an 
inverse "discreteness parameter." As such we should expect that fluctuations 
will be important if (n) is small and the contrary to be true if it is large. We 
shall see that, in fact, this is the case. 

One additional point should be discussed before leaving the question of the 
statistical method. In calculating quantities such as C) /p) and ( ) we face a 
problem with having the statistical quantities in the denominator. While the 
operator ( ) Equation (2), can be applied in principle it would be prohibitively 
difficult to do in practice. The escape fron. this difficulty is effected by defining 
bp = p - (p ) and writing 



We would expect this series, though most likely not a convergent one, to be a 
useful asymptotic one if fluctuations are small, i.e., bp/(r') << 1. As we might 
expect this series turns out to be expressible, after some rearranging, as a 
power series in 1 (n), our "discreteness parameter." 
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ill. Galactic Propagation Model 

We now turn to the question of what particular form the response function 
p(X - X i . t - t .) shall have. The form that this function will take is completely 
dependent on the model one chooses for the propagation of cosmic rays through 
the galaxy. In principle one could choose as complicated a model as one wished 
and perform the integrals indicated in Equation (3) numerically. 

Since, however, we are interested primarily in determining how the overall 
results are affected by certain general parameters of a model we shall consider 
here a model that is relatively simple and amenable to calculation. Further- 
more, it is my belief that this model is not without a physical basis. It is in 
fact suggested by the picture of galactic cosmic ray propagation proposed by 
Parker^, 15, 16 in which the cosmic rays are relatively free to stream along the 
galactic spiral field all the while leaking out of the surface of the galactic disk. 
This leakage is accomplished by means of the "bubble blowing" instability also 
discussed by Parker. 14 

We shall therefore adopt the following picture of cosmic ray propagation in 
the galaxy. When a supernova explodes it suddenly fills a bubble of volume L 0 3 
with a hot cosmic ray gas where L 0 is of the order of the galactic disk thick- 
ness. This bubble is threaded by spiral (on the average) magnetic field of the 
galactic disk and is thus able to expand only along the field lines. (This is why 
we have considered response functions of only one variable.) This bubble 
expands in both directions along the field direction behind a "front" that moves 
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with velocity V while behind the front the density and pressure are uniform. 
While the expansion is taking place the total number of particles is decaying 
exponentially with a time constant r c due to leakage out of the sides. This 
picture can be thought of as propagation down a leaky pipe. From these con- 
siderations we may write the density or response function as 


-(t-t i )/r c 


p(x-x,. t-t,) = r Lo ' v ( t - t rf u(t-t.)u(v(t- tl ) +L 0 - |x-xj) (5) 

where U(X) is the Heaviside step function and K is a normalizing coefficient. 
From considerations of continuity we have 


Vi (X. t) = 


v(x-x t ) 

[ L o + V(t- ti )] 


( 6 ) 


and 


j . - v . p. . 

These formulae hold true for any point in the galaxy that is on a field line 
that treads the initial bubble , for any other point there will be no effect at all. 

There are, of course, many other possibilities for models and response 
functions; diffusive models, or models with more detailed hydrodynamic 
expansions for example. However, I believe that the model described above 
exhibits certain features that would have to be present in any model. The 
minimum length L 0 has no effect on the average density (p) , however, (/O 2 ) and 
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higher order moments diverge as L 0 - 0. This simply means that the proba- 
bility of a supernova exploding arbitrarily near the point of observation has an 
overwhelming effect on all fluctuation phenomena. A diffusion model would have 
to have a built in maximum velocity of propagation otherwise the well known 
infinite propagation speed at t = 0 characteristic of pure diffusion solutions 
causes ( j 2 ) and higher terms to diverge. Finally, I believe that leakage out of 
side of the "pipe" is physically called for by the arguments of Parker. 

IV. Calculations 

I have used Equations (5) and (G) in evaluating* the integrals of Equation (3). 
The series expansion of Equation (4) has been carried out through all terms of 
order 1 , (n). In the calculations the free expansion velocity V has been consid- 
ered a variable parameter and the dependence of (r), (r) + yi(Sr 2 ) , (r) 

- »^>. <v>. (v) * and (v) - y^(&v 2 ) upon this parameter are shown 

in Figures 1 and 2. 

If we consider the RMS deviations to be a measure of the uncertainty of any 
prediction of an observable from a particular model we see that the results for 
the mean lifetime of cosmic ray particles is not too spectacular. (In this model 
lifetime is identical with lifetime in the disk since halo storage is not consid- 
ered.) We see that the best one can do is about 6% accuracy but this is 
insignificant compared to present uncertainties. 

*Data used for the galactic configuration were: Length of spiral flux tube (arm) — 92 Kpc, Earth 
36 Kpc from inner end, L 0 = 0. 1 Kpc, event rate ~ 0.2 per year for our galaxy, galactic radius 
12.5 Kpc, - JO 6 years. 
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However, when we come to the bulk streaming velocity we find a different 
matter entirely. VVc see that the presently observed limit for the streaming 
velo city of £500 km/sec is compatible with a range of inherent streaming 
velocities from one fourth tue speed of light on down. From thi° we can *ee 
that the observations of this streaming velocity may not tell us too much about 
the relevant parameters of a particular propagation model. 
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FIGURE CAPTIONS 


Figure 1. Mean particle age (^) and (r) plus and minus the R.M.S. deviation 
y A dr2 ) as a function of the expansion velocity V. 

Figure 2. Bulk streaming velocity (v) and (v) plus and minus the R.M.S. 


deviation y(c>v 2 ) as a function of the expansion velocity V. 
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